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Abstract 

We present an extensive analysis of AF = 2 observables and of B — > X s j in 
the framework of a specific Maximally Gauged Flavour (MGF) model of Grinstein 
et al. including all relevant contributions, in particular tree- level heavy gauge boson 
exchanges whose effects are studied in detail in the present paper for the first time. 
The model allows in principle for significant deviations from the Standard Model 
predictions for AMB ds , mixing induced CP-asymmetries S^k s and and 
B — > X s ^ decay. The tension between ek and S^Ks present in the SM can be 
removed by enhancing \ek\ without modifying S^k s - In this case, we find that in 
this model i) the results for S^^ and B — > X s j turn out to be SM-like, ii) the exclusive 
determination of \V u b\ is favoured and most importantly hi) the values of AM^ and 
AMb s being strongly correlated in this model with ek turn out to be much larger 
than the data for the central values of input parameters: AMb a ~ 0.75/ps and 
AMb s ~ 27 /ps. Therefore, from the present perspective, the model suffers from 
a serious ejf-AM^ tension. However, this tension can be softened considering 
theoretical and parametric uncertainties and in particular the decrease of the weak 
decay constants. On the other side, the model can be strongly constrained considering 
the theoretically cleaner ratios AMs d / AMb s and BR(B + — > t + v) / AMB d and we 
find that it is unable to remove simultaneously all the SM tensions on the data. 
Finally, we compare the pattern of flavour violation in MGF with selected extensions 
of the SM. 
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1 Introduction 



The Standard Model (SM) of Particle Physics is successful in describing particles and 
their electroweak and strong interactions, still, several aspects are problematic. In this 
paper, we concentrate on the Flavour Problem. 

The introduction of additional symmetries beyond the SM gauge group acting on the 
three fermion generations can produce realistic mass hierarchies and mixing textures. The 
Lagrangian is invariant under the gauge group of the SM and under the additional flavour 
symmetry at an energy scale equal or higher than the electroweak one. Fermion masses 
and mixings arise once these symmetries are broken, spontaneously or explicitly. Such 
flavour models differ from each other in the nature of the symmetries and the symmetry 
breaking mechanism. On the other hand, they all share the same top-down approach: 
the main goal is the explanation of fermion masses and mixings by the introduction of 
flavour symmetries; only as a second step their phenomenological consistency with FCNC 
processes (sometimes) is investigated (see Refs. [1,2] and references therein). 

A bottom-up approach consists in first identifying a low-energy effective scheme in which 
the contributions to FCNC observables are under control and subsequently in constructing 
high-energy models from which the effective description can be derived. The so-called 
Minimal Flavour Violation (MFV) [3-6] follows this second approach. The fact that so 
far no evident deviations from the SM predictions have been found in any flavour process 
observed in the hadronic sector [7] , from rare decays in the kaon and pion sector to B decays 
at supers-factories, can be a sign that any physics beyond the SM does not introduce 
significant new sources of flavour and CP violation with respect to the SM. In Refs. [8-11], 
this criterion has been rigorously defined in terms of flavour symmetries, considering an 
effective operator description within the SM. More in detail, restricted to the quark sector, 
the flavour symmetry coincides with the symmetry of the SM Lagrangian in the limit of 
vanishing Yukawa couplings. This symmetry can be written as the product of non-Abelian 
SU{3) terms, 

G f = SU(3) Ql x SU(3) Ur x SU(3) Dr , (1.1) 

and three additional U(l) factors, that can be arranged to correspond to the Baryon 
number, the Hypercharge and a phase transformation only on the right-handed (RH) down- 
type quarks. Interestingly, only the non-Abelian terms of Gf control the flavour structures 
of the quark mass-matrices, while the U(l) factors can only be responsible for overall 
suppressions [11]. The 5 , C/(2)^-doublet Ql and the 5'[/(2) Z/ -singlets U R and D R transform 
under Gf as 

Q L ~(3,1,1), U R ~ (1,3,1), L>*~(1,1,3). (1.2) 

In order to write the usual SM Yukawa terms, 

C Y = Q L y d D R H + Q L y u U R H + h.c. , (1.3) 

where H = 1T2H* , manifestly invariant under Gf, the Yukawa couplings are promoted to 
dimensionless fields - called spurions - with non-trivial transformation properties under 
G f : 

34 ~ (3, 3,1), y d ~(3,l,3). (1.4) 
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Following the MFV ansatz, quark masses and mixings arise once the electroweak symmetry 
is spontaneously broken by the Higgs VEV, (H) = vj \/2 with v = 246 GeV, and the spurion 
fields obtain the values, 

ft I m d \ /2/ m «00\ 
y d = — m s and y u = — V f m c , (1.5) 
V \ m b J v \ m t ) 

where V is the unitary CKM matrix. 

Recently, several papers [12-14] appeared where a MFV-like ansatz is implemented in 
the context of maximal gauge flavour (MGF) symmetries: in the limit of vanishing Yukawa 
interactions these gauge symmetries are the largest non-Abelian ones allowed by the La- 
grangian of the model. The particle spectrum is enriched by new heavy gauge bosons, 
carrying neither colour nor electric charges, and exotic fermions, to cancel anomalies. Fur- 
thermore, the new exotic fermions give rise to the SM fermion masses through a See-Saw 
mechanism, in a way similar to how the light left-handed (LH) neutrinos obtain masses by 
the heavy RH ones. Moreover, the MFV spurions are promoted to scalar fields - called 
flavons - invariant under the gauge group of the SM, but transforming as bi-fundamental 
representations of the non-Abelian part of the flavour symmetry. Once the flavons de- 
velop suitable VEVs, the SM fermion masses and mixings are correctly described. Still, 
Refs. [12-14] do not provide a natural mechanism for the specific structure of the flavon 
VEVs. This mechanism is highly model dependent, as discussed in Refs. [15,16], in con- 
trast to the fermion and gauge sectors. Such scalar fields may have a phenomenological 
impact, but it is above the scope of the present analysis to provide a realistic explanation 
for the flavon VEV alignment and we will therefore not include these scalar contributions. 

Even if this approach has some similarities to the usual MFV description, the presence of 
flavour-violating neutral gauge bosons and exotic fermions introduces modifications of the 
SM couplings and tends to lead to dangerous contributions to FCNC processes mediated 
by the new heavy particles. Consequently, the MGF framework goes beyond the standard 
MFV and a full phenomenological analysis of this NP scenario is mandatory to judge 
whether it is consistent with all available data. 

In this paper we focus on the specific MGF realisation presented in Ref. [12], even if our 
analysis can be easily applied to other models with gauge flavour symmetries. In particular, 
we extend the study performed in Ref. [12] and point out that the parameter space of such 
a model can be further constrained performing a full analysis on meson oscillations. The 
number of parameters is much smaller than in other popular extensions of the SM and 
therefore it is not obvious that the present tensions on the flavour data can be removed or 
at least softened. Indeed, we observe that the model, while solving the Ek — S^Ks tension, 
cannot simultaneously remove other SM flavour anomalies, which in some cases become 
even more pronounced. 

Relative to Ref. [12] the new aspects of our analysis are: 

- In addition to new box-diagram contributions to AF = 2 processes, considered al- 
ready in Ref. [12], we perform a detailed analysis including the tree-level exchanges 
of new heavy flavour gauge bosons. These diagrams generate LR operators that are 
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strongly enhanced, by the renormalisation group (RG) QCD running, relatively to 
the standard LL operators and could a priori be very important. 

- The impact of the new neutral current-current operators, arising from integrating 
out the heavy flavour gauge bosons, to the B — > X s ^ has been studied in Ref. [17] 
and we apply those results to the model. 

- We point out that for a value of \V u b\ close to its determination from exclusive decays, 
i.e. in the ballpark of 3.5 x 1CT 3 , the model can solve the present tension between Ek 
and S^k s - For slightly larger values of \V u b\, the model can still accommodate the 
considered observables within the errors, but for the inclusive determination of \V u b\ 
it suffers from tensions similar to the SM. 

- We scan over all NP parameters and present a correlated analysis of Ek, the mass 
differences AM Bjs the B + -> t + v and B -)• X s j decays, the ratio AM Bd /AM Bs , 
the mixing-induced CP asymmetries S^k s and S^, and the b semileptonic CP- 
asymmetry A b sl . 

- We find that large corrections to the CP observables in the meson oscillations, Ek, 
S^Ks and S^, are allowed. However, requiring e k to stay inside its 3a error range, 
only small deviations from the SM values of S^k s and are allowed. 

- We find that requiring £k-S^k s tension to be removed in this model implies the values 
of AM^ and AM Ba to be significantly larger than the data. While the inclusion of 
theoretical and parametric uncertainties and in particular the decrease of the weak 
decay constants could soften this problem, it appears from the present perspective 
that the model suffers from a serious Ek — AMb s ,cL tension. 

- We also investigate the correlation among two theoretically cleaner observables, 
AMs d / AMb s and BR(B + — > t + v) / AMs d - In this way, we strongly constrain the 
parameter space of the model and conclude that the tension in BR(B + — > t + u), 
present already within the SM, is even increased. 

- We compare the patterns of flavour violation in this model with those found in the 
original MFV, the MFV with the addition of flavour blind phases and MFV in the 
left-right asymmetric framework. 

- As a by-product of our work we present a rather complete list of Feynman rules 
relevant for processes in the quark sector. 

The structure of the paper is shown in the table of contents. 



2 The Model 

In this section we summarise the relevant features of the MGF construction presented 
in Ref. [12], dealing only with the quark sector. The flavour symmetry is that of eq. (1.1), 
but it is gauged. The spectrum is enriched by the corresponding flavour gauge bosons and 
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by new exotic quarks, necessary to cancel the anomalies: in particular the new quarks are 
two coloured RH S , f/(3)g i -triplets, one LH SU (3) u R -tri-p\et and one LH SU(3) ^-triplet. 
In table 1, we list all the fields present in the theory and their transformation properties 
under the gauge groups. 
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1 


1 


3 
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3 
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SU(3) Dr 
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1 
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Table 1: The transformation properties of the fields under the SM and flavour gauge sym- 
metries. 

With this matter content, the most general renormalisable Lagrangian invariant under 
the SM and flavour gauge groups can be divided into three parts: 

£ = £kin + £-int ~ V [H, Y u , Y d ] . (2.1) 

The first one, Ckin, contains the kinetic terms of all the fields and the couplings of fermions 
and scalar bosons to the gauge bosons. The covariant derivative entering accounts 
for SM gauge boson-fermion interactions and additional flavour interactions involving new 
gauge bosons and fermions: 

D,D J2 *9fN f (A f ),, (A,)„ = £(^) M ^- , (2.2) 

f=Q,U,D a=l 

where gf are the flavour gauge coupling constants, Nf the quantum numbers, AJ the flavour 
gauge bosons and A^^ the Gell-Mann matrices. 

The second term in eq. (2.1), Ci n t, contains the quark mass terms and the quark-scalar 
interactions: 

Ant = A u Q L H ty UR + X' U ^ UL Y U ^ UR + M u ^ UL U R + 

_ _ _ (2.3) 

+ \ d Q L H q dR + X'Jf dL Y d ^ dR + M d ty dL D R + h.c. , 

where M Ujd are universal mass parameters and A^ d are universal coupling constants that 
can be chosen real, through a redefinition of the fields. 

The last term in eq. (2.1), V [H, Y u , Y d ], is the scalar potential of the model, containing 
the SM Higgs and the flavons Y u d . The mechanisms of both electroweak and flavour 
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symmetry breaking arise from the minimisation of this scalar potential. It has not been 
explicitly constructed in Ref. [12] and it is beyond the scope of the present paper to provide 
such a scalar potential (see Ref. [16] for a recent analysis). Therefore, we assume that the 
spontaneous breaking of the electroweak symmetry proceeds as in the SM through the 
Higgs mechanism and that the spontaneous flavour symmetry breaking is driven by the 
flavon fields Y u d which develop the following VEVs: 



(Yd) 



Yd: 



(Y U ) = Y U V. 



(2.4) 



Here Y u ,d are diagonal 3x3 matrices and V is a unitary matrix. We emphasise that, despite 
the similarity to eq. (1.5) of MFV, the matrix V is not the CKM matrix and the vacuum 
expectation values (Yu,d) do not coincide with the SM Yukawa matrices. This is illustrated 
by moving to the fermion-mass eigenbasis. In what follows we focus on the up-quark sector, 
but analogous formulae can also be written for the down-quark sector. The LH and RH 
up-quarks mix separately giving rise to SM up-quarks u l RL and exotic up-quarks u'^^. 




° u (R,L)i S 'U(R,L)i 
■ Su (R,L)i Cu (R,L)i 




(2.5) 



where c U{RL)i and s U(RL)i are cosines and sines, respectively. Denoting with mp the mass of 
the up-type /' = {u l , u n } quark, what follows is a direct inverse proportionality between 
m„i and m„n: 



m u i m u n = M u X u — . 



We can express these masses in terms of the flavour symmetry breaking parameters: 



m,, 



S U Ri S U L i 



-K(Y U ) 



(2.6) 



(2.7) 



URi ULi 

where a straightforward calculation gives 



-"URi 



X u v m u H — y/2 M u m u i 



\/2 M u m u n — X u v m u i 



2 

m a 



m 



(2.8) 

These results are exact and valid for all quark generations. However, taking the limit 
m u H ^> m u i, we find simple formulae that transparently expose the behaviour of the pre- 
vious expressions. In this limit we find 



v X U M U 
^KiYuV 
m„i A« v 



a u Li ~ 



°URi ~ 



U i V2M U 



m 



(2.9) 



(2.10) 



as it is in the usual see-saw scheme in the limit of (Y u )i 3> M u , v. These simplified relations 
are valid for all the fermions, apart from the top-quark for which the condition m# ^> m t 
is not satisfied and large corrections to eq. (2.10) are expected. 
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From eq. (2.10) we see that to reproduce the correct SM quark spectrum, Y u must have 
an inverted hierarchy with respect to the SM Yukawas. 

The presence of new exotic quarks has a relevant impact on the SM couplings. Indeed, 
the charged current-current interactions including SM and heavy quarks are governed by 
a 6 x 6 matrix which is constructed from the unitary 3x3 matrix V of eq. (2.4) and the 
c u Li , Cd Li , s ULi and Sd Li with (i = 1,2,3) introduced above. Adopting a matrix notation, 
the non-unitary 3x3 matrices 

c UL Vc dL , s UL V s dL (2.11) 

describe the charged current-current interactions within the light and heavy systems, 

respectively. The analogous matrices 

Cu L Vs dL , s UL Vc dL (2.12) 

describe the charged current-current interactions between light and heavy fermions. In 
this notation, c u ^ L and s Ut d L are diagonal matrices, whose entries are c Uj d Li and s Ut d Li , 
respectively. Moreover, we point out that in the no-mixing limit, c u ^ Ll —> 1 and s u ^ Li — > 0, 
the (non-unitary) 6x6 matrix reduces to 

(I I) • (213) 

In this case the CKM matrix coincides with the unitary matrix V. As soon as the mixing is 
switched on, the CKM is modified to include c U) d L , which breaks unitarity. However, these 
deviations from unitarity are quite small (see sec. 6.2.2). Moreover no new CP violating 
phases appear in the resulting CKM matrix. At first sight their absence implies no impact 
of new contributions to the CP- violating observables S^k s and S^. However, this is not 
the case due to the modification of the CKM matrix and the presence of flavour gauge 
bosons. In this respect, this framework does differ from the original MFV of Ref. [8]. 

A consequence of the modification of the CKM matrix is the breaking of the GIM 
mechanism if only SM quarks are considered in loop-induced processes. However, once 
also the exotic quarks are included the GIM mechanism is recovered. We return to this 
issue in sec. 3.4. 

The interactions with the Z boson and the Higgs field are modified too. Their effects 
have been already discussed in Ref. [12] and it turned out that the largest constraint comes 
from the modified Zhb coupling. 

Once the flavour symmetry is spontaneously broken by the flavon VEVs, the flavour 
gauge bosons acquire masses and mix among themselves. Using the vector notation for the 
flavour gauge bosons, 

X — (Aq, ■ ■ ■ , Aq, Ay, . . . , Afj, A} Dl . . . , A%) , (2-14) 
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the corresponding mass Lagrangian reads 



and 
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1 

mass 2 


jvi aX 






>«> { A 


( M uu) ab 




"<n> {a 


( M2 nn) ab 






( M Qu) ab 


-( M uo) ba - 2 


( M io) ab 


= ( M DQ) ba = -\ 
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-M 2 , 



DQ 



M DD, 



A SU(3)i Ag[/(3) } 



a \b 

SU(3)i A SU(3) 



a \ 6 

517(3) J A SU(3) 



A SU(3) 



A SU(3) (Y u ) 



A su(3) (Yd) A 6 



SU(3) 



(Yd) 



(2.15) 



(2.16) 



In general, the diagonalisation of this mass-matrix is only numerically possible; for the 
rest of the paper we shall indicate with M.\ the diagonal matrix of the gauge boson mass 
eigenstates A m , where m = 1, ... ,24, and with W(A m , Af), where / = {Q, U, D} and 
a = 1,...,8, the transformation to move from the flavour-basis to the mass-basis (see 
App. A.2). 



3 AF = 2 Transitions 

3.1 Effective Hamiltonian 

In the model in question the effective Hamiltonian for AF = 2 observables with external 
down-type quarks consists at the leading order in weak and flavour-gauge interactions of 
two parts: 

- Box-diagrams with SM IV-boson and up-type quark exchanges. Due to the mixing 
among light and heavy quarks, there are three different types of such diagrams: 
with light quarks only, with heavy quarks only or with both light and heavy quarks 
running in the box, as shown in Fig. 1. If only exchanges of SM quarks are considered, 
the GIM mechanism is broken in these contributions. It is recovered when also the 
exchanges of heavy quarks are taken into account. 

- The tree-level contributions from heavy gauge boson exchanges of Fig. 2, that gen- 
erate new neutral current-current operators, which violate flavour. 

In principle one could consider box-diagrams with flavour- violating neutral heavy boson 
exchanges but they are negligible with respect to the tree-level contributions. 
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u\ v! 1 




u\ u' 1 



Figure 1: The box-diagrams contributing to K° — K° mixing. Similarly for B® — B® mixing. 

d 




— s 




1v,Pl,1h p R 



Figure 2: The tree-diagrams contributing to K° — K° mixing. Similarly, for B® — B® mixing. 
A m is a flavour gauge boson mass eigenstate. 



The effective Hamiltonian for AF = 2 transitions can then be written in a general form 



as 



H 



AF=2 
eff 



(3.1) 



where My/ is the mass of the W^-boson, Qi are the relevant operators for the transitions, 
that we list below, and Cj(/x) their Wilson coefficients evaluated at a scale fi, which will 
be specified in the next section. 

While in the SM only one operator contributes to each AF = 2 transition, i.e. QY LL {M) 
in the list of eq. (3.2), in the model in question there are more dimension-six operators. In 
the absence of flavon exchanges, the relevant operators for the M°-M° (M = K, B^, B s ) 
systems are [18]: 



(3.2) 



QX hh (K) = 


(s Q llM P L d a )^YP L d p ), 


QX h \B q ) 


= {b a ^P L q a ){¥^P L q^) : 


Q^ RR (K) = 


(s a ^P R d a )(s^P R d?), 




= (b a 1 ,P R q a )(¥YP R q P ), 


Q\ K {K) = 


(*V^ a )( s "W), 


Ql R (Bq) 


= (b a ^P L q a )(¥YP R q p ), 


Q™(K) = 


{s a P L d a )(sPp R dP), 


Q\\B q ) 


= (b a P L q a )(¥P R q 13 ). 


where P L , R = 


(lT7s)/2. 







In the next section, we collect the Wilson coefficients of these operators separating the 
contributions from box-diagrams and from the tree-level heavy gauge boson exchanges so 
that 

C^^A^Q + A^Q, (3.3) 



where M — K, B d , B s 
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3.2 Wilson Coefficients from Box-Diagrams 

Keeping in mind the discussion around eqs. (2.11) and (2.12) we introduce the mixing 
parameters: 

A, (JO = V*V id , Xi(B q ) = V* b V ig , (3.4) 

where q = d,s and V is not the CKM matrix but the unitary matrix of eq. (2.4). 

Calculating the usual box-diagrams but including also contributions from heavy fermions 
(see Fig. 1) and corrections to IV-quark vertices according to the Feynman rules in App. A.l 
we find the following contributions to the Wilson coefficients relevant for the K° — K° sys- 
tem at the matching scale \it in the ballpark of the top quark mass 1 : 



where 



Ai(^, K) = (c dL1 c dL2 f X ^ K ) X i( K ) C L c Ilj F ^ x i) ' 

tj=l,2,3 

A 2 ( Mt) k) = ( CdL1 CdL2 f X ^ K ) X >( K ) s L 4, F ^ 4) > 

tj=l,2,3 

A 3 (/i t , K) = (c dm c dL2 f Y X ^ K ) X A K ) C L s l L3 F ( x *> x 'j) + S L c l Lj x i) 

ij=l,2,3 

The arguments of the box-functions F are 

2 



X ■> 



m, 



Mi 



w 



x ■ 



m u/ j 



M 



(3.5) 

(3.6) 
(3.7) 

(3.8) 
(3.9) 



where both i and j run over 1, 2, 3. The loop-function F(xi, Xj) is 



F(xi^ Xj) (4 -\- Xj Xj) I2 (xj, Xj) 8 Xj Xj 1% (xj, x^) 



with 



/^(Xj, Xj^ 



1 



Il\Xii Xj / 



{1-Xi)(l-Xj) 

1 



+ 
+ 



Xj ln(xj) 



x\ ln(^ 



+ (i j) 



+ (i o j 



(3.10) 



(3.11) 



(1-Xi)(l - Xj) 

For the B® — B® mixing we have to replace by B q and Cd i2 by c dL1 c dL3 [c dL2 c dL3 ) in the 
case of q — d (q — s). There are no contributions to other coefficients from box-diagrams. 



1 We explain this choice in the context of QCD corrections below. 
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3.3 Wilson Coefficients from Tree-Diagrams 



Calculating the tree-level diagrams in Fig. 2 with the exchange of neutral gauge boson 
mass-eigenstates A m (m = 1,...,24) we find the following contributions to the Wilson 
coefficient at the high scale which is of the order of the mass of the corresponding 
neutral gauge boson: for the K system we have 



A^cr LL M 



47T 2 * 1 



A (K) C VE« 



4tt 2 



24 

E 



Qt 



Ad 

yR 



ds,m 



-i 2 



An' 



24 

E 



2 g 



ds.m 



ds,m 



ds,m 



(3.12) 
(3.13) 
(3.14) 



where the indices d and s stand for the external quarks d and s, while the index m refers 
to the A m gauge boson mass-eigenstate. The corresponding expressions for the Bd (B s ) 
system are easily derived from the previous ones by substituting ds with db (sb) in the 

indices of the couplings. The explicit expression for the couplings (o d 

App. A.2. 



L,R 



are given m 



3.4 Properties 

We note a few properties: 

- Focussing on eqs. (3.6)-(3.8) and the corresponding expressions in the Bd, s systems, 
for a fixed we have in the box-diagram contributions the combination 

(3-15) 

If all fermion masses were degenerate, this combination would be independent of i,j 
and the unitarity of the matrix V would assure the vanishing of FCNC currents. This 
is precisely what one expects from the GIM mechanism. 

- It is possible to arrange the function J 7 in order to match with the usual notation: 
for the K system we write 

S (x t ) — ► St = (c dL1 c dL2 ) 2 (^33 + - 2-7 r i 3 ) , 
S (x c ) — > SP = (c dL1 c dL2 f (-F22 + Jii - 2JF 12 ) , (3.16) 
S (x c , x t ) — > SiP = {c dL1 c dL2 ) 2 (J*23 + F\i - F13 - Fu) ■ 

For the B q systems we define similar functions S\ q that can be simply derived from 
the previous ones by substituting Cd L1 Cd L2 with Cd L1 Cd L3 {cd L2 Q L3 ) m the case of q = d 
(q = s). In particular the combination of the factors are universal. In order to 
recover the So functions from the expressions it is necessary to take the limit 
in which all the cosines are equal to 1 and all the sines are zero. 
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- The appearance of q and Sj factors introduces in general new flavour dependence, 
implying violation of certain MFV relations even in the absence of new CP-violating 
phases. 

- There are no purely new CP-violating phases in this model, but the CP-odd phase of 
the CKM matrix induces sizeable new effects through new contributions to the mixing 
induced CP-asymmetries S^k s an d S^, in the B® — Bj and the B° s — B° s systems, re- 
spectively. Moreover, similarly to the mass differences AMg d s , new flavour- violating 
contributions affect the parameter ek and there are correlations between the new 
physics contributions to all these observables as we shall see below. 

- The heavy flavour gauge bosons show flavour- violating couplings that can be strongly 
hierarchical: looking at the largest values of these couplings we find 

» (s£«L » (eU) di , (ei,) ct » » (ft*). ■ (3.17) 

An example is presented in App. A. 2 for the lightest gauge boson. This hierarchy is 
due to both the mixings among SM and exotic quarks and the sequential breaking 
of the flavour symmetry encoded in the flavon VEVs, as seen from Eqs. (2.16), (A. 3) 
and (A. 4). 

3.5 QCD Corrections and Hadronic Matrix Elements 

The complete analysis requires the inclusion of the renormalisation group QCD evo- 
lution from the high scales, at which the initial effective Hamiltonians given above are 
constructed, down to low energy scales, at which the hadronic matrix elements are evalu- 
ated by lattice methods. A complication arises in the model in question as several rather 
different high scales are involved, such as the masses of the W-boson Myy, the masses of 
the neutral gauge bosons and the masses of heavy quarks m q n . 

Before accounting for this problem we recall a very efficient method for the inclusion of 
all these QCD effects in the presence of a single high scale, which we denote by //# . Instead 
of evaluating the hadronic matrix elements at the low-energy scale, we can evaluate them at 
Hh, corresponding to the scale at which heavy particles are integrated out. The amplitude 
for M — M mixing (M = K°, B^, B®) at the scale \xn is then simply given by 

f 2 M 2 

A(M — >■ M) = -JL-JZ. J2C?(n H )(M\Q?M\M) , (3.18) 

i,a 

where the sum runs over all the operators listed in eq. (3.2). The matrix element for M — M 
mixing is given by 

(M\Q-(pt H )\M) = InkF&PftM), (3.19) 

where the coefficients P*(M) collect compactly all RG effects from scales below \xn as well 
as hadronic matrix elements obtained by lattice methods at low energy scales. Analytic 
formulae for all these coefficients, P?(B q ) and P?(K), are given in Ref. [19], while the 
corresponding numerical values will be given below for some interesting values of fin- 
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The question then is how to generalise this method to the case at hand which involves 
several rather different high scales. There are three types of contributions for which the 
relevant high energy scales attributed to the coefficients quoted above will differ from each 
other: 

f. The SM box-diagrams involving W^-bosons and the SM quarks. Here the scale is 
chosen to be \i t = 0(m t ). 

2. Tree-level diagrams mediated by neutral heavy gauge bosons, A m . Since we are 
taking into consideration the contributions from all such gauge bosons, we shall take 
as the initial scale for the RG evolution in each case exactly the mass of the involved 
gauge boson. 

3. The only problematic case at first sight are the contributions from box-diagrams 
that involve simultaneously heavy and light particles. Here the correct procedure 
would be to integrate out first the heavy fermions and construct an effective field 
theory not involving them as dynamical degrees of freedom. However, as the only 
relevant contribution comes from the lightest exotic fermion 2 , that is t', whose mass 
is relatively close to mt, we can also here set the matching scale to be fit- As the 
dominant effects from RG evolution, included here, come from scales below Mw, this 
procedure should sufficiently well approximate the exact one. 

Having the initial conditions for Wilson coefficients at a given high scale [in and pro- 
vided also the corresponding hadronic matrix elements at this scale are known, we can 
calculate the relevant M — M amplitude by means of eq. (3.18). As seen in eq. (3.19) these 
matrix elements are directly given in terms of the parameters P^K), P^B^) and P°"(B S ) 
for which explicit expressions in terms of RG QCD factors and the non-perturbative pa- 
rameters B^fii) are given in eqs. (7.28)-(7.34) of Ref. [19]: the /i^ denotes the low energy 
scale and it takes the value 2 GeV (4.6 GeV) for the K system (B q systems). 

The Bf(u. L ) parameters are subject to considerable uncertainties. Exception are the 
B\ LL parameters for which a significant progress has been made in the recent years by 
lattice simulations. In the SM analysis, the RG invariant parameters B\ LL are usually 
considered and denoted by Bk and Bb 9 - We report their values in tab. 3. For completeness 
we recall the values of B\ LL that we extracted from the most recent lattice simulations: 

B\ LL = 0.515(14) , for K system 

B\ LL = 0.825(72) , for B d system (3.20) 

B\ LL = 0.871(39) , for B s system . 

As these parameters are the same for VRR contributions we will combine them together 
with the VLL contributions in the final formula at the end of this section. 

Neglecting the unknown 0(a s ) contributions to Wilson coefficients of the remaining 
operators at the high energy scale, our NLO RG analysis involves only the values of the 

2 This is strictly true only for the B q systems, because in the K system due to the CKM suppressions, 
the contribution from c' may be non-negligible, as accounted for in our numerical analysis. Still, the most 
relevant contribution comes from t' . 
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coefficients P[ R (K), P[ R (B d ) and P^ R (B S ) calculated at fi H . We are not considering the 
intermediate thresholds of exotic quarks, since the smallness of a s and the absence of the 
flavour dependence in the LO anomalous dimensions of the contributing operators render 
the corresponding effects negligible. To obtain these values we need only the values of B[ R 
and £?2 , that we report below in the NDR scheme 3 [20,21]: 

B} R = 0.562(39) (46), B*f = 0.810(41) (31) , for K system 

B\ R = 1.72(4) ) , B^ R = 1.15(3)(±f) , for B d system (3.21) 

B\ R = 1.75(3) (Ig 1 ) , B^ R = 1.16(2)(1|) , for B s system. 

In tab. 2, we show the resulting Pi factors for some relevant values of hh- 



fJ>H 


500 GeV 


ITeV 


3TeV 


lOTeV 


P? ll (»h, K) 
Pf R (fi H , K) 
P? ll (Vh, B d ) 
Pf R (tlH, B d ) 
P? LL (» H , B.) 

Pf R (VH, B 8 ) 


0.392 
-35.7 
0.675 
-2.76 
0.713 
-2.76 


0.384 
-39.3 
0.662 
-2.97 
0.698 
-2.97 


0.373 
-45.0 
0.643 
-3.31 
0.678 
-3.31 


0.363 
-51.4 
0.624 
-3.69 
0.659 
-3.69 



Table 2: Central values of Pi factors for fi H = {0.5, 1, 3, 10} TeV. 

Notice that the LR operators, arising from integrating out the heavy flavour gauge 
bosons, are strongly enhanced by the RG QCD running as can be deduced from the values 
of the Pj factors. A priori, such contributions could be very important. 



3.6 Final Formulae for AF = 2 Observables 

We collect here the formulae we shall use in our numerical analysis. The mixing am- 
plitude M\ 2 (i = K, d, s) is related to the relevant effective Hamiltonian through 

2m K (M«Y = (K°\H A e t 2 \K°) , 2m Bq (M« 2 )* = (^°|^f=X> ( 3 - 2 2) 

with q = d,s. The Ki — K$ mass difference and the CP-violating parameter Ek are then 
given by 

AM K = 2 Re (M£) , e K = jfj,\ Im (M£) , (3.23) 

where ip e = (43.51 ± 0.05)° and tt e = 0.923 ± 0.006 takes into account that (f € ^ it/A 
and includes long distance effect in Imr 12 [22] 4 and ImM 12 [24]. The mixing amplitude 

3 These values can be found in Refs. [20,21], where B{ R (B^ R ) is called B 5 (_B 4 ). 

4 This value has been confirmed by lattice and presented with a smaller error in Ref . [23] . 



14 



entering the previous expressions can be decomposed into two parts, one containing the 
LL and RR contributions and the second only the LR ones: 

M* = {M*) x + {M*) 2 , (3.24) 

where 

(^il)i = ™ K [Bk Vi Xl(K) S™ + Bk V2 \l(K) S<*>+ 



+ 2B K r, 3 X 2 (K) \ 3 (K)S£ ) + 

+ p? LL (fi H , k) (a? cr LL M + a? } c™^*; 



(3.25) 

Analogously, for the the B® s — B® s systems the two parts of the mixing amplitude are 
given by: 



(mi = a \^ LF k m B q [v B B Bq Xl(B q ) 



+ p? ll (» h , B q ) (a^c? ll m + a^c vrr m) 



(Mhh = F l m Bq Pt R ^n, B q ) A™Cf*(^) . 

(3.26) 

Here 771,2,3,5 are known SM QCD corrections given in tab. 3 and P?([ih,M) describe the 
QCD evolution from fi B down to [II for the considered system. For the B q systems, it is 
useful to rearrange the definition of the mixing amplitude M\ 2 as follows [25] 

M? 2 = (M? 2 ) SM C Bq e 2i ^, (3.27) 

where C Bda and ip Bda account for deviations from the SM contributions. Therefore, the 
mass differences turn out to be 

AM Bq = 2\M q 12 \ = (AM Bq ) SM C Bq (q = d,s), (3.28) 

where 

«)sm = I (OsmI^ > ( M i 2 )sm = I MsmI^ • ( 3 -29) 

Here the phases /3 ~ 22° and /3 S — —1° are defined through 

V£ M = \V£ M \e-* and V t s s M = -\V t s s M \e~^ . (3.30) 

The coefficients of sin(AM£ d t) and sin(AMs s t) in the time dependent asymmetries in 
Bd ~~ ^ an d Bg — > if><f) are then given, respectively, by: 

Stpk s = sin(2/3 + 2<p Bd ) , ^0 = sin(2|/3 s | — 2^^) . (3.31) 

Notice that in the presence of non- vanishing <p Bd and (p Bs these two asymmetries do not 
measure /3 and j3 s but ((3 + (p Bd ) and (\(3 S \ — (f Bs ), respectively. 
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3.7 The Ratio AM Bd /AM Bs and the B + -¥ r + u Decay 

The expressions for the mass-differences recovered in the previous section are affected 
by large uncertainties, driven by the decay constants F Bds . To soften the dependence of 
our analysis on these theoretical errors, we consider the ratio among AM^ and AM Ba , 
that we call Ram b , and the ratio among the branching ratio of the B + — > t + v decay and 
AM Bd , that we name Rbr/am- 

Indeed, when considering Ram b , we notice that the SM theoretical errors are encoded 
into the parameter £ = 1.237 ± 0.032, that is much less affected by uncertainties with 
respect to the mass differences. When considering the NP effects, we obtain 

(AMgJjM C Bd (aao\ 

Ram b - (AM v — yr~ ■ K 6 - 6Z ) 

{£AM Ba )SM <^b s 

On the other hand, in the SM, the B + — > t + v decay occurs at the tree-level through 
the exchange of the M^-boson. Therefore, the expression for its branching ratio is only 
slightly modified in our model: 

BR(B+ r + u) = G F mB+m r f 1 _ J^V F l + \ Culi Vub c dL fr B+ , (3.33) 

07T \ m B+ / 

where the NP effects are represented by the cosines. Notice that heavy flavour gauge boson 
contributions could contribute only at the loop-level and can be safely neglected, since they 
compete with a tree-level process. Furthermore, in the previous expression we have assumed 
the SM couplings for the leptons to the VT-boson: even if we are not considering the lepton 
sector in our analysis it is reasonable to assume that any NP modification can be safely 
negligible, as these couplings are strongly constrained by the SM electroweak analysis. 

In the ratio Rbr/am, the dependence on F Bd , which is indeed the main source of the 
theoretical error on AM Bd , is cancelled [26,27]: 

3 7rr B+ cl Li c 2 dL3 m2 \V ub \ 2 f m 2 \ 2 
Rbr/am - ^ — — — — TWl^Tr [2 1 ~ > \ 6 - 6 ^) 



4 t]b B Bd So(xt) C Bd M^\V*V td \ 2 \ m 2 B 

where the second fraction contains all NP contributions and we took m B + m Bd , well 
justified considering the errors in the other quantities. The SM prediction of this observable 
should be compared with the data 

Rbr/am = (3.25 ± 0.67) x 10" 4 ps . (3.35) 

See tab. 4 for the SM prediction. 

4 The b semileptonic CP-asymmetry 

In the B g systems, apart from AM Bq , S^k s an d S^, a third quantity providing infor- 
mation on the meson mixings is the b semileptonic CP-asymmetry A b sl [28,29]: 

A b sl = (0.594 ± 0.022) a d sl + (0.406 ± 0.022) a s sl , (4.1) 
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where 



''si 



SM 



12) SM 



(Mi 



12 )SM 



sin0 d = (5.4 ± 1.0) x 1(T 3 sin. 



sin0 s = (5.0 ± 1.1) x 10~ d sin 



(4.2) 



with 



arg(-(M^) 5M /(r? 2 ) 5M 
arg(-(M 1 s 2 ) 5M /(r i2 ) 

SM 



-4.3° ± 1.4° 



0.22° ±0.06 c 



(4.3) 



In the presence of NP, these expressions are modified. Since we have already discussed 
the NP effects on Mf 2 in the previous sections, we focus now only on T q l2 - It is useful to 
adopt a notation for T\ 2 similar to the one in eq. (3.27) for M\ 2 : 



1 12 



(r?2)sM C Bq e 



(4.4) 



where C B is a real parameter. With such a notation we get, 



SM 



(Mf 2 ) 



SM 



c. 



B a 



c. 



Ba 



sin (<p d + 2ip Bq + 2(p Bq ) 



(4.5) 



Notice, that in the MGF context we are considering, the phase (p Bq is vanishing, while C Bq 
is mainly given by c^ L2 Cd Lb Cd Lq ~ 1. As a result the only NP modifications are provided 
by the NP contributions on M\ 2 . 



5 The B X s 7 Decay 

5.1 Effective Hamiltonian 

The decay B — > X s ^ is mediated by the photonic dipole operators Q77 and Q' 7 ^ and 
through mixing also by the gluonic dipole operators Q 8 q and Q' 8G - In our conventions they 
read 

g 

Qt-y = 17^? m b s» P R b a , 

6 (5.1) 
Qsg = ^ m ^a ^ p R T ^ h G% 

and the corresponding primed dipole operators are obtained by substituting Pr with P^. 

The effective Hamiltonian for h — > 37 at a scale /i in the SM normalisation and consid- 
ering only the dipole operators reads 



4G F 

"7^ 



v; s v tb 



AC 7 ^)Q 7l + AC 8G (/i)Q 8G + 

+ Ac* 77 0*)g / 77 + AC 8O ( A *)g , 8O ' 



(5.2) 
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We have kept the contributions of the primed dipole operators Q' 7y and Q' 8G even though 
their Wilson coefficients are suppressed by m s /m), with respect to the unprimed Wilson 
coefficients. However, the mixing of neutral current-current operators into Q' 7l and Q' 8G 
can affect ACy 7 as shown in Ref. [17]. 

Similarly to the Hamiltonian for the AF = 2 transitions, the Wilson coefficients in the 
Hamiltonian can be separated into two parts: 

- The SM-like contribution from diagrams with jy-bosons with modified couplings to 
both SM and exotic quarks of charge +2/3, denoted below by u and u', respectively: 



u , u 




1,9 




The contribution of heavy neutral gauge bosons exchanges with virtual SM and exotic 
quarks of charge —1/3, denoted below by d and d', respectively: 




1,9 



The first contribution has already been considered in Ref. [12], while the second, the 
impact of the heavy neutral gauge bosons on b — > 57, has been recently pointed out in 
Ref. [17]. In particular it has been found that the QCD renormalisation group effects in 
the neutral gauge boson contributions can strongly affect the branching ratio of B — > X s , -f 
and cannot be neglected a priori. 



5.2 Contributions of VF-exchanges 

For the H^-exchange the matching is performed at the EW scale, fiw The Wilson 
coefficients are the sum of t and t' contribution, since d and v! contributions are suppressed 
by their small couplings to b and s quarks. Hence, the Wilson coefficients of Q 7l and Q$g 
are 

A W CV 7 (M = c dL2 c dL3 (cl L3 C^ M (x t ) + sl LS C™{x> t )) , (5.3) 
A w C 8G (nw) = c dL2 c dL3 {c 2 UL3 C^(x t ) + s 2 UL3 C^(x' t )) , (5.4) 



with 



^'"°|^ h - Tir (5 ' 5) 

Q/*i2 ,-vt 3 ,7-1 2 O™ 
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being the SM Inami-Lim functions [30]. 

At last we need to evolve A w C([i\y) down to fif, to obtain the contribution of W 
exchanges to the branching ratio of B — >■ X s 7. The QCD analysis, that involves the SM 
charged current-current operators Q\ and Q2 as well as the QCD-penguins Q3 to Qe, which 
mix with and Qsg below fi t , is the same as in the SM and we proceed as in Ref. [17]. 

5.3 Contributions of A m -exchanges 

The contribution of a neutral gauge boson to the effective Hamiltonian in eq. (5.2) 
derives from integrating out the mass-eigenstate of the heavy flavour gauge boson A m 
at its mass-scale fijj. The Wilson coefficients A^C^/i^) and AaC^q(ijh) have been 
calculated within a generic framework in Ref. [17]. The results for the special MGF case 
we are discussing are fixed by the couplings of the flavour gauge bosons to both SM and 
exotic fermions. Applying general formulae of Ref. [17] to the present case we found that 
these contributions are below 1% and can be safely neglected. As discussed in Ref. [17] 
the reason for such suppression are the See-saw-like couplings of flavour gauge bosons to 
both SM and exotic fermions and the heavy neutral gauge boson masses. 

6 Numerical Analysis 

Having at hand the analytic expressions derived in the previous sections, we are ready 
to perform a numerical analysis of the MGF model in question. 

The first question we ask is whether the model is able to remove various anomalies 
in the flavour data hinting the presence of NP. Since the number of parameters is much 
smaller than in other popular extensions of the SM like SUSY models, LHT model, RS- 
scenario and models with left-right symmetry, it is indeed not obvious that these anomalies 
can be removed or at least softened. We briefly review the flavour anomalies as seen from 
the SM point of view. 

6.1 Anomalies in the Flavour Data 
6.1.1 The e K — S^ Ks Anomaly 

It has been pointed out in Refs. [22,31-33] that the SM prediction for ek implied by the 
measured value of S^k s = sin 2/3, the ratio Ram b , and the value of \V c b\ is too small to agree 
well with the experiment. We obtain the SM Ek value by taking the experimental value 
of S^k S i the ultimate |V^| determination, the most recent value of the non-perturbative 
parameter Bk [34-39], and by including long-distance effects in ImFi2 [22] and ImMi2 [24] 
as well as recently calculated NNLO QCD corrections to Ek [40,41]. We find 5 \ek\ = 
(1.82 ± 0.28) x 10 -3 , visibly below the experimental value. 

On the other hand sin 2/3 = 0.85 ± 0.05 from SM fits of the Unitarity Triangle is 
significantly larger than the experimental value. This discrepancy is to some extent caused 
by the desire to fit both e k [22,31,32] and BR(B + t+u) [42]. 

5 The small discrepancy with respect to the value of Ref. [41] comes solely from updated input values. 
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As demonstrated in [22,32], whether the NP is required in e K or S^k s depends on the 
values of 7, \V u b\ and \V c b\. The phase 7 should be measured precisely by LHCb in the 
coming years while \V u b\ and \V c b\ should be precisely determined by Belle II and Super- B 
provided that also the hadronic uncertainties will be under a better control. 

6.1.2 The |T4 6 |-Problem 

There is a tension between inclusive and exclusive determinations of \V u b\. This means 
that if we take the unitarity of the CKM matrix as granted and also consider the good 
agreement of the ratio Ram b with the data, the inclusive and exclusive determinations 
imply different patterns of NP in CP- violating observables. Indeed one is lead to consider 
two limiting scenarios: 

Scenario 1: Small |V^b|. Here \V u b\ is in principle the exclusive determination, 

\V ub \ = (3.38 ± 0.36) x 10~ 3 . (6.1) 

Within the SM, when the AM^/AM^ constraint is taken into account, one finds 
S^k s ~ 0.67 in agreement with the data, but Ek ~ 1-8 x 10~ 3 visibly below the data. 
As discussed in Refs. [22,32], a sizeable constructive NP contribution to Ek would 
not require an increased value of sin 2/3 relative to the experimental value of S^k s - 
NP of this type would then remove the Ek — S^k s anomaly in the presence of the 
exclusive value of \V u b\. 

Scenario 2: Large |Vub|. In this case \V u b\ corresponds to its inclusive determination, 

\V ub \ = (4.27 ± 0.38) x 10~ 3 . (6.2) 

In this scenario the SM predicts Ek ~ 2.2 x 10~ 3 , in agreement with the data, while 
S^Ks ~ 0-81 is significantly above the data. As discussed in Refs. [22,31], a negative 
NP phase (ps d in B® — B® mixing would solve the Ek — S^k s anomaly in this case (see 
eq. (3.31)), provided such a phase is phenomenologically allowed by other constraints. 
With a negative (pB d , sin 2/3 is larger than S^k s , implying a higher value on \er\, in 
reasonable agreement with data and a better Unitary Triangle fit. 

In both scenarios, new physics contributions to other observables, such as AMs da , are 
expected and a dedicate analysis is necessary. In fact as we will see below the correlations 
between e k and AM Bd s are powerful tests of the ability of MGF to describe properly all 
data on AF = 2 observables. 

6.2 Input Parameters and the Parameter Space of the Model 

Before proceeding with our numerical analysis, it is necessary to fix the input parameters 
and to define the parameter space of the model. 
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G F = 1.16637(1) x 10- 
M w = 80.399(23) GeV 
sin 2 W = 0.23116(13) 
a(M z ) = 1/127.9 
a s {M z ) = 0.1184(7) 



GeV 



[43] 
[43] 
[43] 
[43] 



m u (2 GeV) = 1.7 -r 3.1 MeV [43] 
m d (2 GeV) = 4. 1 -=- 5.7 MeV [43] 



m s (2GeV) = 100^{jMeV 



[43] 



m c {m c ) = (1.279 ± 0.013) GeV [44] 
m 6 (m 6 ) = 4.19±S:51 GeV [43] 
M t = 172.9 ±0.6 ±0.9 GeV [43] 



m K = 497.614(24) MeV [43] 

F K = 156.0(11) MeV [36] 

B K = 0.737(20) [36] 

n e = 0.94(2) [24] 

cp e = (43.51 ± 0.05)° [22] 

% = 1.87(76) [41] 

m = 0.5765(65) [45] 

77s = 0.496(47) [40] 



[43] 
[43] 
[36] 
[36] 
[36] 
[36] 

[36] 

[36] 
[36] 
[45,46] 

r B± = (1641 ±8) x 10~ 3 ps [43] 



= 5279.5(3) MeV 
: 5366.3(6) MeV 
205(12) MeV 
250(12) MeV 
1.26(11) 
1.33(6) 



m Bd = 
m Bs = 
Fs d = 

Fb ° = 

B B d = 

B Bs = 

F Bd ^B Bd = 233(U) MeV 

F Bs y/B~ Bs = 288(15) MeV 
f = 1.237(32) 
Vb = 0.55(1) 



\V US \ ■ 

\Vcb\ ~- 
I T/incl. 
I V ub 
I T /excl 
I V ub 

7 



0.2252(9) [43] 
(40.6 ± 1.3) x 10~ 3 [43] 
= (4.27 ±0.38) x 10~ 3 [43] 
= (3.38 ±0.36) x 10" 3 [43] 
(73t|)° [43] 



Table 3: Values of experimental and theoretical quantities used throughout our numerical 
analysis. Notice that mj(mj) are the masses rrii at the scale rrii in the MS scheme. M t is 
the pole top-quark mass. 



6.2.1 Input Parameters 

In Table 3 we list the nominal values of the input parameters that we will use for the 
numerical analysis, except when otherwise stated. At this stage it is important to recall 
the theoretical and the experimental uncertainties on some relevant parameters and on the 
observables we shall study. 

Considering the K° — K° mixing, remarkable improvements have been made in the case 
of the CP-violating parameter Ek, where the decay constant Fk is known within 1% accu- 
racy. Moreover the parameter Bk is known within 3% accuracy from lattice calculations 
with dynamical fermions [34] and an improved estimate of long distance contributions to 
Ek reduced this uncertainty down to 2% [22,24]. The NNLO QCD corrections to 771 and 
773 [40,41] allowed to access the remaining scale uncertainties that amount according to [41] 
to roughly 6%, dominantly due to the uncertainty in 771. Including also parametric uncer- 
tainties, dominated by the value of \V c b\, Brod and Gorbahn estimate conservatively the 
present error in Ek to amount to roughly 15% [41]. The reduction of this total error down 
to 7% in the coming years appears to be realistic. Further reduction will require progress 
both in the evaluation of long distance contributions and in rji. AM^ is very accurately 
measured, but is subject to poorly known long distance contributions. 

Regarding the — B® mixings, lattice calculations considerably improved in recent 
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years reducing the uncertainties in F Bs 6 and F Bd and also in y B Bg F Bs and y B Bd F Bd 
down to 5%. This implies an uncertainty of 10% in AM Bd and AM Bg within the SM. 
On the other hand, the mixing induced CP-asymmetries and S^k s have much smaller 
hadronic uncertainties. 

The hadronic uncertainties in the ratio Ram b are roughly at the 3% level; the theoretical 
error on the b semileptonic CP- asymmetry A b sl is around the 20% level; the theoretical 
uncertainties in the rate of the B — > X s ^ decay are below 10%; given the 5% uncertainty 
on F Bq decay functions, the branching ratio for B + — > t + u has a theoretical error around 
the 10% and a large parametric error due to \V u b\. 

We stress that the situation with other Bi parameters, describing the hadronic matrix 
elements of AF = 2 operators absent in the SM, is much worse. Here a significant progress 
is desired. 



On the experimental side, Ek, AM Bd , AM Bs and the ratio Ram b are very precisely 
measured with errors below the 1% level. S^k s is known with an uncertainty of ±3% 
and the rate for the branching ratio of the B — > X s 7 decay is known within 10%. On 
the contrary, larger experimental uncertainties affect the measurements of S^^ in DO and 
LHCb, which differ from one another by an order of magnitude, but are still in agreement 
within the 2<r-level, due to the large errors of the single determinations. A sl has only been 
measured by DO its experimental error is around the 20%. Similarly, BR{B + — > r + v) is 
plagued by the same uncertainty. 



6.2.2 The CKM Matrix 

To evaluate the observables we need to specify the values of the CKM elements. As 
already stated in sec. 2, the CKM matrix in this model is not unitary and is defined by 

V = c UL Vc dL , (6.3) 

where V is by construction a unitary 3x3 matrix and ci u ,d) L are the cosines encoding 
the mixing between SM and exotic fermions. From eqs. (2.8) and (2.10), we deduce that 
c (u,d) L ~ 1; except for t and t'. As a result, within an excellent accuracy the CKM matrix 
reads 

(V u d V us V u b \ 
Vcd V cs V cb . (6.4) 
c UL3 Vtd c UL3 V ts c UL3 V t b J 

In this approximation, the deviation from the unitarity of the CKM matrix is 

{V ] V) . . = Sa - 4 3 V* V tj , (V ft) . . = 5ij - s 2 UL3 5 lt Sjt , (6.5) 

The deviations are present only when the top-quark entries are considered and are propor- 
tional to s UL3 . All other entries of the CKM matrix coincide with the corresponding entries 
of the unitary matrix V up to negligible corrections. 



6 Recently a remarkably precise value for Fb b was reported in Ref. [47]: Fb s = (225 ± 4) MeV. Still, 
we shall adopt a conservative approach and use the value of Fb b = (250 ± 12) MeV in our analysis except 
when explicitly stated. 
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The important implication of the latter finding is that the angle 7 in the unitary triangle 
is unaffected by such deviations. In the approximation of eq. (6.4), 



and thus 7 does not depend on c UL3 or s UL3 . 

We state now how we fix the values of the CKM elements. From the tree-level exper- 
imental determinations of \V US \, \V c b\, \V u b\ and 7, we fix the corresponding parameters of 
V. In this way, also the corresponding parameters of the V matrix are univocally fixed and 
using the unitarity of V, we evaluate all the other entries of V. With all entries of V fixed 
we compute the masses and mixings of all fermions and flavour gauge bosons by means of 
eqs. (2.6)-(2.16). Finally, knowing c ULS , we also determine the elements of the third row 
of V. 

6.2.3 The Parameter Space of the Model 

Having determined V what remains is the calculation of the spectrum and the couplings 
of NP particles. In principle they are fixed once, in addition to the SM parameters, we 
fix the seven NP couplings d , Qq , gu , go and the two mass parameters M u and M d in 
eqs. (2.2) and (2.3). Still, their actual determination is subtle since the energy scale at 
which the see-saw relations of eqs. (2.6) hold is a priori not known. We identify this scale 
with the mass of the lightest flavour gauge boson. 

We fix the spectrum and the see-saw scale iteratively using the condition that all exotic 
masses are above m t (m t ). As a first step we evaluate the see-saw relation at m t (m t ) to 
obtain a rough estimate of the masses of exotic fermions and lightest gauge boson. With 
this initial spectrum we run the masses of the SM fermions to the newly defined see-saw 
scale including all intermediate exotic fermion thresholds. The evaluation of the see-saw 
relation corrects the NP spectrum. We repeat the procedure until the values of exotic 
fermion masses and see-saw scale no longer change. Lastly, we evolve the exotic fermion 
masses down to the EW scale. 

For the numerical analysis it is necessary to scan the parameter space of the model. 
We choose \ Ut d G (0, 1.5] and all other couplings {X' ud , gQ, gu, go} G (0, 1.1] to stay 
in the perturbative regime of the theory. The two mass parameters are varied between 
M u G [ 100 GeV , 1 TeV ] and M d G [30 GeV , 250 GeV ] following the discussion in Ref. [12] . 
Unphysical points of the parameter space, namely cases with s u ^ d or c Uyd larger than 1, are 
not considered. Larger M u and M d values decouple the NP from the SM and are therefore 
phenomeno logically irrelevant. With respect to the analysis of Ref. [12] we are scanning 
over all NP parameters, including X' ud and gQ, gu, gu,. 



To present the features of the MGF model we are discussing, we use \V US \ and \V c b\ at 
their central values in tab. 3 and 




(6.6) 



6.3 Results 



3.38 x 10 



3 



and 



7 = 68' 



o 



(6.7) 
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which are among the favoured values within the SM when the experimental values of both 
S^k s and Ram b are taken into account. 

With this CKM matrix we list in tab. 4 the central values for the SM predictions of the 
observables under consideration together with their experimental determinations. 



SM predictions for exclusive 


v ub 




Experimental values 


AM Bd = 0.592 ps" 1 
AM Bs = 20.28 ps^ 1 
Ramb = 2.92 x 10- 2 
SfK s = 0.671 
= 0.0354 

AM K = 0.4627 x lO^ps" 1 
\e K \ = 1.791 x 10- 3 
A b sl = -0.0233 x lO- 2 
BR(b — >■ 57) = 3.15 x 10~ 4 
BR(B + t+v) = 0.849 x 10" 4 
Rbr/am = IA3x 10" 4 ps 


AM Bd = 0.507(4) ps" 1 [43] 
AM Bs = 17.77(12) ps -1 [43] 
Ram b = (2-85 ± 0.03) x 10- 2 [43] 
S fKs = 0.673(23) [43] 
0f = 0.55^;1 [48,49] 
= 0.03 ±0.16 ±0.07 [50] 
AM K = 0.5292(9) x lO^ps" 1 [43] 
\e K \ = 2.228(11) x 10~ 3 [43] 
A b sl = (-0.787 ± 0.172 ± 0.093) x 10~ 2 [28] 
BR(b 57) = (3.55 ± 0.24 ± 0.09) x 10^ 4 [43] 
BR(B + ->■ t + v) = (1.65 ± 0.34) x 10~ 4 [43] 
Rbr/am = (3.25 ± 0.67) x 10~ 4 ps 



Table 4: The SM predictions for the observables we shall consider using the exclusive 
determination of \V u b\ and the corresponding experimental values. 



Comparing these results with the data we make the following observations 

- \ek\ is smaller than its experimental determination, while S^k s is very close to the 
central experimental value, as it should be for the chosen \V u b\ and 7. 

- The mass differences AM Bd and AM Bs are visibly above the data; also their ratio 
Ram b 1S above the experimental determination, but in agreement at the 3<r level. 

- BR(B + — > t + u) is well below the data and consequently also the ratio R B r/am turns 
out to be below the measured central value by more than a factor of two. Even if the 
experimental error in BR(B + — > t + v) is large, the parameter space of the model is 
strongly constrained. Furthermore, from the correlation among Ram b and Rbr/am 
it is evident that the model can only deteriorate the SM tension in these observables. 

- Concerning S^, the predicted value is consistent with the most recent data from 
CDF, DO and LHCb. 

- A b sl is well below the DO data. 

- Finally the predicted central value for BR(B — > Xg'-y) is smaller than the central 
experimental value but consistent with it within the 2a error range. 

Any NP model that aims to remove or soften the anomalies listed above should simul- 
taneously: 

1. Enhance \ek\ by roughly 20% without affecting significantly S^k s - 
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2. Suppress AM Bd and AM Ba by roughly 15% and 10%, respectively. 

3. Slightly suppress Ram b by 3%. 

4. Strongly enhance Rbr/am by 130%. 

5. Moderately enhance the value BR(B — > Xgj) by 5 — 10%. 

As we shall see below, the model naturally satisfies requirements 1., 3. and 5. On the 
other side, it fails in 2. and 4.: indeed the mass differences AMs ds can only be enhanced 
with respect to the corresponding would-be SM value; this enhancement is predicted to 
be significant if ones requires to solve the \sk\-S^k s anomaly. Furthermore, the predicted 
value for the Rbr/am can only be decreased, resulting in a tension on this observables 
more serious than in the SM. 

In what follows we will look closer at the pattern of flavour violations in the MGF still 
keeping the input parameters at their central values. Subsequently we will comment on how 
some of our statements are softened when hadronic uncertainties in the input parameters 
are taken into account. 

Considering now the NP contributions within MGF, we find the following pattern of 
effects: 

- If we neglect the contributions of flavour gauge bosons and of exotic quarks, i.e. 
considering only the would-be SM contributions, the mixing amplitudes M{ 2 and 
BR(B — > X s -f) are reduced with respect to the SM ones due to the modification 
of the CKM matrix, encoded in the mixings c u ^ L - As a result, once the third-row 
entries of the CKM matrix are involved, the would-be SM values of the considered 
observables are smaller than the values reported in tab. 4. 

- The RR flavour gauge boson contributions are negligible for all observables in all the 
parameter space. We shall not consider such contributions in the following descrip- 
tion. 

- \ek\ is uniquely enhanced by the new box-diagram contributions involving exotic 
quarks, while it is uniquely suppressed by heavy gauge flavour boson contributions. 
Among the latter, the LR contributions are the dominant ones, while the LL ones 
are safely negligible. 

- AMB ds are also uniquely enhanced by the new box-diagram contributions, but are 
mostly unaffected by heavy flavour gauge boson contributions. This is in particular 
true for AM^, while for AMg s the latter contributions can be non-negligible either 
enhancing or suppressing it. This is best appreciated when considering the ratio 
Ram b - this observable does not show any dependence on the new box-diagram con- 
tributions, since in MGF the operator structure in box-diagram contributions does 
not change with respect to the SM and the NP effects are the same in the Bd and 
B s systems. As a result any NP effect in this ratio should be attributed to the heavy 
gauge flavour boson contributions, both LL and LR. 
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- The mixing induced CP-asymmetries S^k s an d are unaffected by the new box- 
diagram contributions. Similarly to Ram b > this allows to see transparently the heavy 
gauge flavour boson contributions, which was much harder in the case of AM# ds . 
We find that S$k s * s om Y affected by LL contributions and can only be suppressed. 

depends on both LL and LR contributions. Interestingly, the NP contributions 
interfere destructively with the SM contribution such that the sign of S,^^ can in 
principle be reversed in this model. Similar conclusions hold for A b sl : it is not affected 
by box-diagram contributions, the LR contributions are almost completely negligible 
and the LL ones are the only relevant enhancing \A b sl \ towards the central value of 
the experimental determination. 

- Finally, the branching ratio of B — > Xg'j can be significantly affected by the modifi- 
cations in the SM magnetic penguin contributions that can only enhance this observ- 
able, as already pointed out in Ref. [12]. The heavy gauge flavour boson contributions 
are negligible as discussed in Ref. [17]. 

Having listed the basic characteristic of NP contributions in this model we will now present 
our numerical results in more detail stressing the important role of correlations among 
various observables identified in this model by us for the first time. 

6.3.1 Correlations Among the Observables 

In this section we discuss correlations among the observables. They will allow us to 
constrain the parameter space of the model and see whether this model is able or not to 
soften, or even solve, the anomalies in the flavour data. 
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Figure 3: The correlation of Ek and S^Ks- The shaded grey regions are the experimental 
la -3a error ranges, while the cross is the central SM values reported in tab. 4- The colour 
of the points represent the percentage of the box-diagram contributions (purple) and of the 
flavour gauge boson ones (red) in Ek- In the NP points the theoretical error on Ek is 
included. 
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In fig. 3a, we show the correlation between Ek and S^k s - The plot confirms that the 
exclusive value of \V u b\ is favoured in this model. Indeed the NP contributions are able 
to solve the £k—S$k s anomaly in a reasonably large region of the parameters space. This 
happens when the Ek prediction approaches the data due to box-diagram contributions 
(purple points), while S^k s is mostly unaffected. This is in particular possible when the 
flavour gauge bosons contributions are negligible. When the flavour gauge boson contribu- 
tions are significant (red points) S$k s is uniquely suppressed relatively to the SM value. 
However, as seen in the figure a combination of large box contributions as well as flavour 
gauge boson contributions (yellow points) allows bringing ek in agreement with the data 
while keeping S^k s within the 2a experimental error range. 

On the other hand, points for which the box contributions are negligible and instead 
the flavour gauge boson contributions dominate in e k (purely red points) cannot explain 
the observed value of Ek- However, this kind of contributions are best suited for the case 
with the inclusive determination of \V u b\, reported in fig. 3b. Still, there exist no points, 
which simultaneously bring ek and S^k s in a 3a agreement. The flavour gauge bosons 
contributions are not large enough to suitably correct S^k s - We conclude that the inclusive 
determination of \ V u b\ is disfavoured in this model and we shall not further pursue this case. 

In fig. 4, we present the correlations between Ek and AMs ds - From (a) and (b), we 
conclude that the model cannot solve the \sk\ ~ S^k s anomaly present in the SM, without 
worsening the already moderate agreement of this model with the AM^ s experimental 
data. Indeed the values AMs d ss 0.75/ps and AMb s ~ 27/ps are so larger that in the case 
that the central values of the weak decay constants Fg . s do not change in the future, but 
their and other input parameter uncertainties are further reduced, we will have to conclude 
that the model fails to describe the AF = 2 data. 

However, we should emphasise that this problem could be avoided if the values for the 
weak decay constants Fs da are smaller than the ones used in the plots (a) and (b). Indeed, 
in (c) we adopt a 15% reduced value for F Bd , close to its 3a value, while in (d) the last 
determination of Fb s reported in Ref. [47]. This input, modifies the SM values to be 

AM Bd = 0.43 ps" 1 and AM Bs = 16.4 ps" 1 , (6.8) 

such that now the enhancements of these observables by NP is welcomed by the data. From 
plots (c) and (d) we deduce that the NP contributions and the requirement of agreement 
of \ek\ with data within MGF automatically enhance AM^ s . Even if also in this case 
AMs d s are found above the data, the model is performing much better than in the previous 
case. This exercise shows that on one hand it is crucial to get a better control over hadronic 
parameters in order to obtain a clearer picture of NP contributions and on the other hand 
that other more precise observables should be analysed until the uncertainties on AMs d s 
are lowered. 

In fig. 5, we show the correlation between the ratio of the mass differences, Ram b , and 
the ratio among the branching ratio of the B + — > t + u decay and the AM Bd , Rbb,/am- 
Both observables have negligible theoretical uncertainties and are therefore very useful to 
provide strong constraints on the parameter space. On the left, the red points correspond 
to an ek prediction in agreement with the data at 3a, while the blue ones do not satisfy 
the Ek constraint. This plot largely constrains the parameter space of the model; for only 
very few points Ek, Ram b and Rbr/am agree at 3cr-level with the data simultaneously. 
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(c) F Bd = 175 MeV (d) F Bs = 225 MeV 



Figure 4: Correlation plot of Ek with AM^ and AM Bg on the left and right, respectively. 
In the upper plots, we use for F Bd s the values reported in tab. 3, while for the lower plots 
we adopt smaller values: F Bd is reduced down to 85% of its value, close to the 3a error 
level and F Bg is taken to be the new determination reported in Ref. [47]. In the plots above 
we have not included the la error in ek to best illustrate the interplay of box- (purple) and 
tree- contributions (red). 

Furthermore, all red points correspond to values for Rbr/am smaller than the SM 
prediction and therefore the model can only worsen the SM tension. In the case that the 
experimental sensitivity to the BR(B + — > t + v) improves, it will be possible to further 
constrain and possible exclude the present MGF model. 

Having reduced the parameter space, we concentrate now on a few other predictions of 
the model. In fig. 6, we show the mt> — Mj^a parameter space, where mt> is the mass of 
the exotic partner of the top-quark and Mj$a the mass of the lightest neutral gauge boson; 
the corresponding particles have the best chances to be detected at the LHC. The red and 
blue points are those identified in fig. 5 to agree and disagree in Ram b , Rbr/am and Ek 
with the data at the la-3a level, respectively. 

Interestingly, the phenomenological results presented above hold not only for light but 
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Figure 5: Correlation plot for R&m b ana " Rbr/am- The grey regions refer to the experi- 
mental la-3a and 2a-3a error ranges for Ram b and Rbr/am, respectively. The big black 
point refers to the SM values reported in tab. 4- On the left, red (blue) points refer to 
agreement (disagreement) of the points prediction of Ek and the data at 3a level. On the 
right the colours represent contribution of boxes and trees in Ek ■ 
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Figure 6: m t / — M424 parameter space. For all red points in the plot Ram b , Rbr/am an d 
Ek agree with the data at 2a -level. 

also for heavy i"s. Also the mass of the lightest flavour gauge boson is not bounded. 

Furthermore, in fig. 7 we show two correlation plots which represent also clear predic- 
tions for this model. Plot (a) is the correlation for and A b sl showing that only tiny 
deviations from the SM values are allowed: this turns out to be an interesting result for 
S^, which is indeed close to the recent determination of LHCb. On the other hand, A b sl 
has only been measured by DO, but hopefully LHCb will also have something to say in the 
near future. Once the experimental uncertainties are lowered, such clear predictions will 
be essential to provide the final answer on how well this model performs. 

In plot (b), we show the correlation of BR(B — > X s j) and m' t . We confirm the finding 
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Figure 7: Correlation plot of S^^ and A b sl on the left and BR(B — > X s 7) and m' t on the 
right. Grey regions refer to the experimental error ranges. The big black point refers to 
the SM values reported in tab. 4- In r ^d the points for which Ram b , Rbr/am and Ek agree 
with the data at 3<r level, in blue all others for which there is no agreement. 

of Grinstein et al. that in this model the NP contributions to BR(B — > X s / y) always 
enhance it towards the central experiment value. However, interestingly only very small 
enhancements of this branching ratio are allowed when also the bounds from AF = 2 
observables are taken into account. 



7 Comparison with other Models 

A complete comparison of the patterns of flavour violation in MGF with corresponding 
patterns found in numerous models [51] would require the study of AF = 1 processes, 
however already AF = 2 observables allow a clear distinction between the MGF and the 
simplest extensions of the SM . Here we just quote a few examples: 

- In the original MFV framework restricted to LL operators, the so-called constrained 
MFV [52], the \ek\ —S^Ks anomaly can only be solved by enhancing \ek\ since in this 
framework S^k s remains SM-like. In this framework then only the exclusive value of 
\V u b\ is viable. An example of such a framework is the model with a single universal 
extra dimension (UED) for which a very detailed analysis of AF = 2 observables 
has been performed in [53]. In fact this is a general property of CMFV models as 
demonstrated in [54]. Thus after \ek\ has been taken into account and contributions 
from tree-level heavy gauge boson exchanges have been eliminated MGF resembles 
CMFV if only AF = 2 processes are considered. However AF = 1 processes can 
provide a distinction. In fact whereas in MGF the NP contributions uniquely enhance 
BR(B — > Xs'j), in UED they uniquely suppress this branching ratio [55]. Concerning 
the |ex|-AMs ds tension MGF and CMFV are again similar. 
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- The 2HDM framework with MFV and flavour blind phases, the so-called 2HDMp^py 
[56], can on the other hand be easily distinguished from MGF. In this model NP 
contributions to Ek are tiny and the inclusive value of \V u b\ is required in order to 
obtain the correct value of \sk\- The interplay of the CKM phase with the flavour 
blind phases in Yukawa couplings and Higgs potential suppress S^k s simultaneously 
enhancing the asymmetry S^. As in the case of MGF this asymmetry is SM-like or 
has a reversed sign. It is together with the value of |V^&| which will distinguish 
MGF from 2HDM MW . 

- Finally, let us mention the left-right asymmetric model (LRAM) for which a very 
detailed FCNC analysis has been recently presented in [57]. As this model has many 
free parameters both values of \V U b\., inclusive and exclusive, are valid. The model 
contains many new phases and the \sk\ — S^k s anomaly can be solved in many 
ways. Moreover, the model struggles with the Ek constraint due to huge neutral 
Higgs tree-level contributions. However, as demonstrated in Section 7 of that paper 
a simple structure of the right-handed mixing matrix gives a transparent solution to 
the \ek\ — S^k s anomaly by enhancing \er\, keeping S^k s at the SM value and in 
contrast to MGF automatically suppressing AMs ds and significantly enhancing S^^. 
While MGF falls back in this comparison, one should emphasise than on the MGF 
has very few parameters and provides the explanation of quark masses and mixings, 
while this is not the case in the LRAM. 



8 Conclusion 

We have presented an extensive analysis of AF = 2 observables and B — > A s 7 for a 
specific MGF model presented in Ref. [12], which is of special interest due to the small 
number of new parameters. In particular we performed a detailed study of the effects of 
tree-level contributions due to the presence of heavy flavour gauge bosons. 

Our main findings are as follows. The model predicts a clear pattern of deviations from 
the SM: 

- Enhancements of \ek\ and AM^ ds in a correlated manner by new box-diagram con- 
tributions and suppression of \ek\ by tree-level heavy gauge boson contributions with 
only small impact on AM^ s . 

- Mixing induced CP-asymmetries S^k s and are unaffected by box-diagram con- 
tributions, but receive sizeable destructive contributions from tree-level heavy gauge 
boson exchanges such that the sign of can be reversed. However, these effects 
are basically eliminated once the Ek constraint is taken into account. 

- The ek — S^k s anomaly present in the SM is removed through the enhancement of 
\ek\, leaving S^k s practically unmodified. This is achieved with the help of box- 
diagram contributions in the regions of the parameter space for which they are domi- 
nant over heavy flavour gauge boson contributions, which interfere destructively with 
the SM amplitudes. 
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- This structure automatically implies that in this model the exclusive determination 
of \V u b\ is favoured. 

- b semileptonic CP-asymmetry A sl , that a priori could receive large contributions 
from the tree-level flavour gauge boson diagrams, remains close to the SM values 
once requiring to be in agreement with the data. 

- Most importantly, the Ek constraint implies the central values of AMg d s to be 
roughly 50% higher than the very precise data. This disagreement cannot be cured 
fully by hadronic uncertainties although significant reduction in the values of Fb a a 
could soften this problem. 

- We have pointed out that the ratio of the mass differences, AM^/AMb s and the 
ratio of the B + — > t + v branching ratio and AMs d , together with ek, provide strong 
constraints on the parameter space of the model. Furthermore, the correlation among 
these two observables encodes a serious tension on the flavour data that can only be 
deteriorated in the model. 

- In agreement with Ref. [12], we find that BR(B — > X s 7) is naturally enhanced in 
this model, bringing the theory closer to the data, still only small corrections are 
allowed by the AF = 2 bounds. 

- We have demonstrated how this model can be distinguished by means of the flavour 
data from other extensions of the SM. 

In summary, the great virtue of this model is its predictivity, such that within the 
coming years it will be evident whether it can be considered as a valid description of low- 
energy data. Possibly the most transparent viability tests of the model are the future values 
of Fe d s and \V u b\. For the model to accommodate the data, Fs d s have to be reduced by 15% 
and \V u b\ has to be close to 3.4 x 10 -3 . Violation of any of these requirements will put this 
extension of the SM in trouble. We emphasise that the triple correlation between AM^ s , 
Ek and S^k s was instrumental in reaching this conclusion. The same result is derived 
by the complementary study on the Rbr/am — Ram b correlation. Our work shows that 
the study of correlations among flavour observables and the accuracy of non-perturbative 
parameters like B^ and Fs d3 are crucial for the indirect searches for physics beyond the 
Standard Model. 
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A Feynman Rules for MGF 

A.l Couplings of SM Gauge and Goldstone bosons 



7 coupling 

The photon coupling remains unchanged; 
proportional to the quark charges Q u and 



G coupling 

The gluon coupling remains unchanged; 
proportional to the colour generators A^^. 



W ± coupling 
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H coupling 
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with the actual values for /, / and C: 
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• \ 

1A u S u RiS u Li 


Godidi 


Cr = 


-C l 


= +i\dSdRiCdLi 


Godi 


Cr = 


-C l 


= —iXdCdRiSdLi 


God^d; 


Cr = 


-C l 


= —iXdCdRiCdLi 


God'^i 


Cr = 


-C l 


= +i\dSdRiSdLi 



G ± coupling 



G ± 




= i(P L C L + PrCr) 



f 



with the actual values for /, / and Cl- 



G + Uidj 


Cr 




+ ^dC u LiVijSdRj 




C l 




SuRi Vij CdLj 


G+u 1 ^ 


Cr 




— ^dSuLiVijCdRj 




C l 




+ ^uC u Ri Vij SdLj 


G+Uid'j 


Cr 




— ^dCuLiVijCdRj 




C l 




+ ^uCuRiVijCdLj 


G+u>d 3 


Cr 




+ ^dSuLiVijSdRj 




C l 




SuRi Vij $dLj 


G~djUi 


Cr 




— ^u s uRiV*jCdLj 




C l 




+ ^dC u LiVi*jSdRj 




Cr 




+^uCuRiV*jSdLj 




C l 




— ^dS u LiV*jCdRj 


G'dju'i 


Cr 




+ ^uC u RiV*jCdLj 




C l 




— ^dCuLiV*jCdRj 


G~d'jUi 


Cr 




— ^uS u RiV*jSdLj 




C l 




+ ^dS u LiV i *SdRj 



A. 2 Couplings of Flavour Gauge Bosons 

There are three types of flavour gauge bosons, Aq, Ajj, and A a Dl which are flavour 
eigenstates but not mass eigenstates. We denote the mass eigenstates with A m , where 



34 



m = 1,...,24. 



8 \T 



X — (^Q; • • • j ^Q; ^[/; • • • j A/> • • • J ^Z?) 



A 24 ) T . 



Flavour basis \ and mass basis <p are connected through the transformation 

X = W<f, (A.l) 
where W is obtained numerically by diagonalising the mass matrix in eq. (2.15) such that: 



Ml = W T MlW 



(A.2) 



and M.a is a diagonal mass-matrix. 
We define: 

U T = (u,c,t,u',c',t') and D T = (d, s,b, d' , s' ,b') 

such that the coupling of the flavour gauge bosons to the quarks are described by the 
Lagrangian-part 

where m is understood to run from 1 to 24 and the tensors Q = Q[Cl,r, 9q,u,d] can be 
read off from the couplings of flavour eigenstates Aq, Ajj, and Ad to the quarks, that are 
listed below: i.e. 



(^l)i3,1 _ ~2 C «ii (^S[/(3))l3 c u L3 > (Or) 42,18 ~~ Sd ^ (^SC7(3))l2 c d i2 ■ ( A - 3 ) 

The rotation to the mass-eigenstates of the heavy gauge bosons redefines the couplings: 



C d + Ua^Gt + GidijjkUj ■ 



where 



&*) ijm = Z>fa m ' X*) fe)^ , with a = M. 



(A.4) 



Notice that in the rest of the paper we either use {i, j} G {1, 
the (SM or exotic) quark flavour: i.e. 



Ql 



L,R 



15, m 



Ql 



L,R) . > 

uc'm 



(sU) 



23,m 



6} or refer directly to 
(A.5) 



'L,R 



sb,m 



A Q coupling 



with the actual values for /, / and C: 



A" 
A Q 



i^(P L C L + PrCr) 



f 
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AgUiUj 


C R = 




C L = 




C R = 




c L = 


Q 1 J 


C R = 




c L = 




C R = 




c L = 


d'ldj 


Cr = 




c L = 




Cr = 


c L = 


J\.Qdidj 


Cr = 


c L = 


A a Q d\dj 


Cr = 




c L = 



Ajj coupling 



+s U Ri(y \ a su ^v^)ijS U Rj 

-\-C u Li(V\g U ^ 3 - ) V^)ijC u Lj 
+CuRi(VXg u ^V^) ij C uRj 

~ s uRi{V X%jj($)V^\ijC u Rj 
+C u Li(VXg u ^V )ijS u Lj 

—c U Ri(VXg u ^V)ijS U Rj 

"I" s uLi ( V^SU(3) ^ ) V C uLj 



+SdRi\X a SU ($y)ijSdRj 
+CdLi{X a su ^)ijC dL j 

~\~ c dRi (-^5(7(3) ) *J C dRj 
+ SdLi{X a su ^)ijSdLj 

— SdRi (tfsufiyij C dRj 
+CdLi{X a su ^)ijSdLj 

— CdRi (^5{J(3) ) ij s dRj 
+SdLi{X a su ^)ijCdLj 




^^(PlCl + PrCr) 



with the actual values for /, / and C: 



AjjUiUj 



Afjiiiu'j 



Cr = ~\~C u Ri{\^jj^)ijC u Rj 

Cr = +St t iji(Agj^/g\)ij - S w iy 
Cl = +C u Li(Xg U ^)ijC u Lj 
Cr = ^C u Ri(\gjj{^)ijS u Rj 
C l SuLi{Xgjj(£^)ijC u Lj 
Cr = ^S u Ri{\ a gjjf^)ijC u Rj 
Cl — —C u Li(X su ^)ijS u Lj 



Ae> coupling 




i^(P L C L + PrCr) 



with the actual values for /, / and C: 



Apdidj 


Cr 


= +CdRi(Xg u{3) 


)ijC d Rj 




c L 


= + s rfLi(^SC/(3) 


^ijSdLj 


w 3 


Cr 


= + S dRi(Xg [J ^ 


) ij SdRj 




c L 


= + C dLi{X a su{3 y 


ij CdLj 


j\.£)didj 


Cr 


= + c dRi(Xg u ^ 


)ijSdRj 




c L 


= —SdLiiXgu^ 


)ijCdLj 


Apd^dj 


Cr 


= + S dRi(Xg u ^ 


)ij c dRj 




c L 


= ~ C dLi(Xgjj^ 


'ij^dLj 



B Couplings of the Lightest Flavour Gauge Boson 

As an example, we report in fig. 8 and 9 the couplings of the lightest gauge boson 
to the light down- and up-type quarks for the exclusive V u b case and scattering on all 
the parameters of the model. The colour-coding corresponds to the mass splitting of the 
lightest to the next-to-lightest flavour gauge boson mass. Namely, red points correspond 
to large splitting while black points to degeneracy, yellow to intermediate splitting. We 
observe that the flavour-violating couplings with the last (first) two generations are in 
general the largest (smallest) ones. This is related to the sequential breaking of the flavour 
symmetry encoded into hierarchical structure of the flavon VEVs. 
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Figure 8: The couplings of the lightest flavour gauge boson to the down-type SM quarks. 
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Figure 9: The couplings of the lightest flavour gauge boson to the up-type SM quarks. 
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